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The rise of the C IV mass density at z < 2.5* 



o 

(N 
£ 

m 

(N 



o 

u 

43 

Oh 
6 

(N 
> 

(N 

d 

On 
O 



13 



Valentina D'Odorico 1 !, Francesco Calura 1 , Stefano Cristiani 1 ' 2 , Matteo Viel 1 ' 2 

1 INAF-OATS, Via Tiepolo 11, 34U3 Trieste, Italy 

2 INFN/National Institute of Nuclear Physics, via Valeric 2, 34-127 Trieste, Italy 



ABSTRACT 

The cosmic evolution of the metal content of the intergalactic medium puts stringent 
constraints on the properties of galactic outflows and on the nature of UV background. 
In this paper, we present a new measure of the redshift evolution of the mass density of 
C iv, ficrVj hi the interval 1.5 < z < 4 based on a sample of more than 1500 C iv lines 
with column densities 10 12 < iV(Civ) < 10 15 cm~ 2 . This sample more than doubles 
the absorption redshift path covered in the range z < 2.5 by previous samples. The 
result shows a significant increase of ficiv towards the lower redshifts at variance with 
the previously pictured constant behaviour. 
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1 INTRODUCTION 

The cosmological mass density of Civ, ficiv, observed as a 
function of redshift is a fundamental quantity closely related 
to the metal enrichment of the intergalactic medium (IGM). 
Its app arent lack of evolution in the redshift interval z ~ 
[1.5, 5] (|Songailall20fjll : |Pettini et al.ll2003l ; lBoksenberg et all 
2003) is puzzling since both the physical conditions of the 
IGM and the properties of the ionizing background are 
thought to evolve between these epochs. 

Remarkable efforts have been spent in recent years 
to extend the measure of f2civ to re dshift larger than 5 
(|Rvan- Weber et al.ir200d ; ISimcoell2006l ) where a decrease of 
the s tar formation rate density is observed |Bunker et all 
2006). If flciv is dominated by the metals produced in situ 
by the observed star forming galaxies, we would expect a 
decrease of its value at those redshifts. Vice versa, the value 
of fiery could remain constant if it reflects the metallicity 
of a diffuse medium pre-enriched at very high redshift. It 
should be noted, however, that this is a simplified scenario 
since, as redshift increases, the observed Civ absorptions 
likely trace gas in structures of decreasing over-density and 
also the ionizing spectrum evolves in shape and intensity. 
As a consequence, the behaviour of ficiv could be different 
from that of VLq a nd of the mean IGM metallicity (see e.g. 
ISchave et af]|2003l ). 

The most recent measurements of C iv absorptions in 
spectra of QSOs at z ~ 6 seem to i ndicate a downturn 
in the C iv mass density at z > 5 (|Becker et all 120091 : 
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iRvan- Weber et al]|2009l ). though based only on 3 detected 
C iv lines. 

At redshift z < 4.5, a f undament a l mea surement of 
ficiv has been carried out by ISongailal (|200ll . SOI). How- 
ever, the redshift interval 1.5 < z < 2 is poorly sampled by 
the considered QSO spectra. A more uniform redshift co ver- 
age is provided by the sample o f lBoksenberg et all (|2003h al- 
though with fewer QSO spectra. Both analysis are consistent 
with a constant behaviour of ficrv in the ran ge [1.5, 4.5]. At 
z < 1 , recent results based on HST UV data (|Cooksev et all 
l200Sh give ficiv = (6±1) x 10" 8 corresponding to a 2.8±0.5 
increase over the 1.5 < z < 5 values. 

In this paper, we present a new measurement of Qciv in 
the redshift range [1.5, 4] based on a sample of 25 high reso- 
lution, high signal-to-noise QSO spectra plus an additional 
sample of 8 QSO spectra from the literature. 

The rest of the paper is organized as follows. The data 
are presented in § 2. In § 3 the analysis is carried out with 
the computation of ficiv- The results are discussed in § 4. 
Throughout this paper, we assume S7 m = 0.26, Qa = 0.74 
and h = i7 /(100km s^Mpc" 1 ) = 0.72. 



2 OBSERVATIONAL DATA SAMPLE 

The core of our sample is formed by the high resolu- 
ti on, high signal-to-n oise QSO spectra al r eady described 
in lSaitta et all l|2007l ) and lD'Odorico et all (|2008l ). Most of 
them were obtained wit h the Ultraviolet an d Visual Echelle 
Spectrograph (UVES) (|Dekker et all l2000h at the Kueyen 
unit of the ESO VLT (Cerro Paranal, Chile) in the frame- 
work of the ESO Lar ge Programme (LP): "The Cosmic Evo- 
lution of the IGM" (|Bergeron et al.ll2004l ). 

In this work, 3 more QSOs were added to that sam- 
ple, mainly to increase the redshift extension above z ~ 3. 
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UVES spectra of the QSO: Q0055-269, PKS2000-330 and 
PKS1937-101, were downloaded from the ESO Archive and 
reduced with the UVES pipeline following the standard pro- 
cedure. The continuum level was determined by interpolat- 
ing with a cubic spline the region of the spectrum free from 
evident absorption features. 

For all the QSO in the sample (see Table [Q), the Civ 
forest was defined as the interval between the Lya emis- 
sion and 5000 km s _1 from the Civ emission to avoid the 
proximity region where most of the intrinsic systems are 
found. The absorption features present in this wavelength 
interval were identified inspecting the spectra by eye to look 
for the most common doublets (Civ, Mgn and Siiv). Then, 
other lines were identified testing their compatibility with 
the Civ, Mgn and Siiv redshifts. Finally, lines whose iden- 
tity was still unknown after this operation, were associated 
with metal systems detected in the Lya forest or recognized 
as part of other multiplets (e.g., Fell). 

The C iv doublets were fitted with Voigt profiles us- 
ing the LYMAN context o f the MIDAS reduction package 
(|Fontana fc Ballestedll995h . A minimum number of compo- 
nents was adopted to fit the velocity profile in order to reach 
a normalized \ 2 ~ 1- The fit parameters for all the detected 
C IV lines are reported in Table [J]. In the following, we will 
refer to C IV components or simply C IV lines meaning the 
velocity components in which every absorption profile has 
been decomposed. However, in order to compare our results 
with previous works and with data at lower resolution, we 
will work also with C IV systems formed by groups of com- 
ponents. C IV systems were defined in the following way: for 
each list of C IV components corresponding to a single QSO 
the velocity separations among all the lines have been com- 
puted and sorted in ascending order. If the smallest sepa- 
ration is less than efo m i n = 50 km s _1 (corresponding to the 
velocity separation adopted by SOI) the two lines are merged 
into a new line with column density equal to the sum of the 
column densities, and redshift equal to the average of the 
redshifts weighted with the column densities of the compo- 
nents. The velocity separation are then computed again and 
the procedure is iterated till the smallest separation becomes 
larger than dvmin- 

Our sample consists of 1023 Civ velocity components 
with column densities 10 12 < JV(Civ) < 10 15 cm -2 and of 
508 C IV systems in the same column density range. 

The C IV absorptions in the spectra of 19 of the QSOs 
forming our samp l e wer e already identified and fitted by 
IScannapieco et all |2006l ) using the software package VP- 
FIT0. We refer to that paper for a careful analysis of the 
clustering properties of C IV lines and a comparison with 
the properties of Siiv, Mgn and Fell lines. Unfortunately, 
the authors did not publish the individual Voigt parameters 
of their fitting. 

2.1 Additional data sample 

In order to further increase the number of C IV lines and to 
extend the sample to higher redshift, we have considered the 
C IV lines fitted in 9 QSO spectra observed with the High 

1 The complete tables are available only in electronic format. 

2 http : / / www.ast.cam. ac.uk/~rfc/vpfit. html 



Table 1. Relevant properties of the QSOs forming the total sam- 
ple. See text for further details. 
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Resolution Echelle Spectrometer (HIRES) at Keck at a res- 
olution and signal-to -noise ratio similar to tho se of our spec- 
tra and reported in iBoksenberg et all l|2003l . BSR03). The 
fit with Voigt profiles was carried out by the authors with 
VPFIT. The main difference between LYMAN and VPFIT is 
that the number of components fitted to a given velocity 
profile is, i n general, larger u sing the latter (see also the dis- 
cussion m ISaitta et alll2007l ). This is seen also in the present 
case, in particular from the comparison of the C IV lines de- 
tected in the spectrum of the QSO B1422+231, which is the 
only object in common between the two samples. We find 
that in all cases the number of components found with VP- 
FIT is larger or equal to that found with LYMAN. However, 
when the total column density of each absorption system is 
considered the difference between the two fitting procedures 
becomes negligible. 

This is shown in Fig. [1] where the number of C IV com- 
ponents (upper panel) and the total column density of C IV 
systems (lower panel) obtained with the two fitting packages 
are compared. While the number of components is signifi- 
cantly larger in the fit by BSR03 the total column densities 
are most of the time in very good agreement. 

The sample of BSR03 is formed by 577 C IV components 
and 302 C IV systems in the column density range 10 12 < 
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Table 2. Civ absorption lines: HE1341-1020 (z c 



2.142) 



1.699561 
1.699690 
1.699818 
1.700051 
1.700284 
1.700996 
1.700812 
1.701957 
1.703648 
1.854894 
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1.910986 
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1.914971 
1.915286 
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2.041588 
2.084978 



b 

km/s 
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29.6 ±2.5 
6.9 ±0.3 
9.5 ±0.4 
10.2 ± 1.0 

7.5 ±0.5 
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10.5 ±0.3 



13.44 
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12.87 
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13.04 
13.00 
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12.33 
11.88 
12.16 
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12.20 
13.06 
12.80 
12.88 
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± 0.06 
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N(Cw) < 10 15 cm -2 . In the following, we will refer to our 
sample of C IV lines as Sample A and to the BSR03 sample 
as Sample B (excluding the lines of B1422+231, which are 
already in Sample A). The total sample is the sum of Sample 
A and Sample B. All the QSOs forming this sample are 
reported in Table Q] with their emission redshift and the 
redshift range covered by the C IV forest. 



3 DATA ANALYSIS 

3.1 The Civ column density distribution 
function, f(N) 

f{N) is extremely sensitive to the velocity decomposition 
of absorption features, since it is defined as the number of 
lines per unit column densi ty and per unit redshift absorp- 
tion path, dX (|Tvtlerll 1987m ." In the assumed cosmology, the 
redshift absorption path is: 

dX = (1 + z) 2 [tl m (l + zf + Q A r 1/2 dz. (1) 

In order to compare with previous works, f(N) has been 
computed for the Civ systems in samples A and B in the 
common redshift interval, 1.6 < z < 3.6 (see Fig. [2]). The 
excess of low column density systems in BSR03 is due to 
their over-decomposition of the C IV velocity profiles with 
respect to our fit. A maximum likelihood fit to the data 
with column densities log iV(Crv) 13 (binned in Fig.[2]for 
display purposes only) to a power law of the form f(N) = 
B N~ a , gives an index a = 1.71 ±0.07 for sample A and a — 
1.8±0.1 for sample B. Both are in agreement with the result 
by S01 based on C IV systems defined in th e same way. These 
are steeper than the lEllison et al.l (|2000r i fit of 1.44 ± 0.05 
based on the very high SNR spectrum of B1422±231 and 
extending the power law down to log iV(Civ) ~ 12.3. 

In order to estimate the incompleteness of our data for 
log TV (C iv) ^ 13 we have performed the following simu- 



Figure 1. Upper panel: Number of components of the Civ ab- 
sorption systems in the spectrum of B1422±231 as fitted by us 
(solid line) and by BSR03 (dashed line). Lower panel: Column 
densities of the C IV systems of B1422±231 computed by us (open 
dots) and by BSR03 (solid triangles). 



Table 3. Results of incompleteness tests 



log N(C iv) Fraction of detected C IV 


12.00 


0.60 


12.30 


0.82 


12.60 


0.97 


12.78 


0.97 


13.00 


1.00 



lations. Four fake C iv doublets have been generated for 
each spectrum in our sample, with column densities 12 ^ 
TV(Civ) ^ 13 and redshift chosen randomly from the red- 
shift ranges reported in Table[T] The b- values of the fake C IV 
lines were drawn at random from the observed distribution 
of b- values in our data. These fake Civ doublets were then 
added to the real spectra and searched for by the member of 
our team who had previously identified the real Civ lines. 
Given the visual, rather than automatic, character of our 
searches, only two such trials were performed. The results 
of these tests are collected in Table [3] 

We found that we could recover essentially all C iv dou- 
blets as long as log N(C iv) ^ 12.6. For log N(C iv) = 12.3 a 
correction factor of 1.2 has been determined (31/38 fake C IV 
systems detected in the incompleteness tests). In the lowest 
column density bin, logiV(Civ) = 12, 24 over 40 Civ lines 
have been identified, corresponding to a correction factor of 
1.7. The average Doppler parameter of the undetected Civ 
lines is (6) ~ 17.5 km s -1 , while detected lines have (b) ~ 9 
km s~ . 
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Figure 2. Comparison of the column density distribution func- 
tion of C IV systems for Sample A (open dots) and Sample B 
(crosses) in the common redshift range 1.6 < z < 3.6. The bin- 
size is 10 ' 3 N(C iv) cm" 2 and the error bars are ±1 a based on 
the number of points in each bin. The dashed and dotted lines 
are power laws of the form f(N) = BN~ a with index a = 1.53 
and 1.71, respectively (see text). 



Extending the range of fitted data to our completeness 
limit, log-/V(Clv) ^ 12.6, the maximum likelihood indexes 
become a — 1 .53 ± 0.04 for sample A and a = 1.61 ± 0.07 
for sa mple B, in agreement with the result bv lEllison et al.l 
(2000). If, the correction factors reported above are applied 
to the points of Sample A, their values increase slightly to 
coincide with the values in Sample B shown in Fig. [2] The 
corrected value of the lowest column density bin is still too 
low to be fitted with the previously derived power laws. 

3.2 The number density of Civ lines, dn/dX 

The evolution with redshift of dn/dX computed for Civ 
lines in Sample A is reported in the upper panel of Fig. [3] 
In particular, we have split our sample into weak and strong 
lines characterized by column densities 12 ^ log N(C iv) ^ 
13 and 13 < log N(C iv) ^ 15, respectively. The two pop- 
ulations show different behaviours, with the weak lines be- 
ing consistent with no evolution and the strong ones signif- 
icantly decreasing in number towa rds high e r reds hifts. The 
latter trend confirms the result by Steidcl (1990). Further- 
more, the ratio of the number of strong and weak lines shows 
a steady decrease with redshift. 

3.3 The two-point correlation function of Civ 
lines, £ (v, v + Av) 

We computed £(v,v + Av) for the weak and strong Civ 
components and for two redshift intervals 1.5 < z < 2.3 



redshift 

Figure 3. Upper panel: Redshift evolution of the number density 
of the C IV lines in Sample A divided into low column density 
(open dots, 12 ^ log N (Civ) ^ 13) and high column density 
(solid dots, 13 < log N (Civ) ^ 15) absorptions. Lower panel: 
Redshift evolution of the ratio between the number of high and 
low column density C IV lines. 



and 2.3 < z < 3.5 (see Fig. [4j. The obtained two-point 
correlation functions show a characteristic clustering scale 
ro (in redshift space) for the sample at high column den- 
sities which is about 1.5 — 2 times larger than for the low 
column density lines (~ 850 vs. 500 km s _1 in the high 
redshift bin and 700 vs. 314 km s _1 for the low redshift 
bin, with la errors of about 150 km s _1 ). Converting these 
scales into comoving Mpc h -1 at a given redshift, it ap- 
pears that the low column density sample has clustering 
properties simi lar to those of Lyman-break galaxies (LBGs) 
at z ~ 3 (e.g. iPorciani fc Giavaliscoll200"2l ). while large col- 
umn density C iv systems are more strongly clustered and 
possi bly sample denser enviro nment and more massive ob- 
jects. lAdelberger et all (|2005l ) reached similar conclusions 
from the analysis of a large sample of LBGs observed in the 
fields of 23 high redshift QSOs. 

3.4 The redshift evolution of the mass density of 

Civ 

We used our sample of absorption lines also to compute 
the mass density of C iv as a fraction of the critical density 
today: 

fi CIV = E^EL f Nf(N)dN, (2) 

Cpcrit J 

where Ho = 100 h km s _1 Mpc~ 1 is the Hubble constant, 
mciv is the mass of a Civ ion, c is the speed of light, 
Pcrit = 1.88 x 10~ 29 /i 2 g cm" 3 and f(N) is the Column 
Density Distribution Function (CDDF). Since f(N) cannot 
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Figure 4. Two-point correlation function of the C IV lines in 
Sample A divided into low column density (open dots, 12 ^ 
logA r (Clv) ^ 13) and high column density (solid dots, 13 < 
log N (Civ) 15) absorptions. Solid lines trace the result for the 
redshift bin 1.5 ^ z ^ 2.3 while dashed lines refer to the bin 
2.3 < z ^ 3.5. 



be recovered correctly for all column densities due to in- 
completeness and poor statistics, the integral in eq. [2] can 
be approximated by a sum: 



Homciv E.A^Civ) 



C p crit AX 

with an associated fractional variance: 
/5ficiv\ 2 _ EJ^(Civ)] 2 
V Qciv ) 



(3) 



(4) 



as proposed by IStorrie-Lombardi et a 3 l|l996l) . Error bars 
have been computed also using a bootstrap technique to 
build 1000 samples of QSO spectra, based on the observed 
sample, for each redshift bin and determining the standard 
deviation of the resulting distribution of S7civ values. This 
estimate is always slightly larger (by a factor of ~ 1.5 or 
smaller) than the one based on eq. [4] and has been used as 
the reference one. 

The value of f2civ significantly depends on the column 
density range over which the sum or the integration are car- 
ried out, and as a consequence on the resolution and SNR of 
the available spectra. To take this aspect into account, we 
have computed three sets of values to be compared consis- 
tently with different data in the literature. 

ficiv obtained for the C IV components in Sample A 
with column densities 12 ^ log (Civ) ^ 15 is compared 
with the values reported in the original paper by BSR03 in 
Fig. [5] Our results confirm and strengthen the increasing 
trend with decreasing redshift of f2civ which was already 
present in the data by BSR03. 



Figure 5. Comparison of the SlciV determinations based on Civ 
lines in Sample A (open dots) and the results by BSR03 (crosses). 
Error bars are 1 a. 
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Figure 6. f2civ determined for the systems in the total sample 
(open dots) with 1 a error bars from bootstr ap method compared 
with the results by S01 (solid squares) and bv lPettini et al.l d2003T ) 
(empty triangle). 
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Table 4. f2civ Ior the systems in the total sample selected in the 
reported column density intervals 



10 



> i 

u 

~5 



3 A 
Redshift. 



Figure 7. Estimates of SlciV m the redshift range z ~ 
— 6 for Civ syste ms with 13.8 < logJV(Clv) < 15: cross 
ICooksev et"ai] d2009h : open dots total sample this work; s olid tri- 
anales iPettini et al j <2003i) and iRvan- Weber et al.l j2009T) ; 95 % 
confidence interval lTjecker et al.l J2009t) ! 



In order to carry out the comparison with the data 
set by SOI, the most used in the literature for the red- 
shift range [1.5,4.5], Slciv has been computed from the 
C IV system^j of the total sample with column densities 
12 < IogJV(Crv) ^ 15. In Fig. [6] these results are reported 
together with the value of SOI corrected for the assumed 
cosmology (they assumed Einstein-de Sitter) and with error 
bars derived from their plot and transformed to la. The 
present data sample represents an increase in the absorp- 
tion redshift path of a factor 2.5 and 2 in the redshift bins 
[1.5,2.0] and [2.0,2.5] respectively, with respect to S01. In- 
deed, our estimate in the lowest redshift bin is ~ 4.5 a larger 
than the previously accepted value, suggesting a trend of in- 
creasing Civ mass going to lower redshifts. The computed 
values with the associated errors are reported in Tabled 

We have not applied corrections to the values of ficiv 
due to the incompleteness of our observations for C IV lines 
with column densities < 12.6, since we estimated them to 
be less than 3 % (less than 1 % in the lowest redshift bin). 
This is due to the small correction factors determined in Sec- 
tion 3.1 and to the small contribution of low column density 
lines to the final values. 

Recently, the determination s of Qciv have bee n ex- 
tended to very low (z < 1, e.g. Cooksev et alj 200S ) and 
very high (z > 5 lRvan- Weber et al.l2009l ; lBecker et alj 



2009) 



redshift. Observations in these redshift ranges are more dif- 



3 We remind the readers that C IV systems have been built adopt- 
ing a minimum velocity separation of 50 km s _1 as reported in 
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ficult since the C IV transition moves to the UV and to the 
IR region of the electromagnetic spectrum, respectively. As 
a consequence, the estimates of ficiv are based on lower 
resolution and lower SNR spectra limiting the detectability 
of C IV lines to larger column densities. In order to study 
the evolution of ficiv in the whole redshift range between 
z ~ and 6, we have carried out a third computation of 
the values of f2civ for C IV systems in the column density 
range 13.8 ^ log AT (Civ) ^ 15, changing the definition of 
the redshift bins in order to have comparable absorption red- 
shift paths covered in each one of them. The column density 
lower boundary is the limit to which the higher (and lower) 
redshift surveys cited above are sensitive (being based on 
lower resolution and SNR spectra than those analysed in 
the present work). The results are shown in Fig. [7] where 
the po int at z ~ 4.69 is the determination bv lPettini et al.l 
(2003) corrected for the con sidered column density range 
by IRvan- Weber et all (|2009T l. The plot shows clearly the 
smooth growth of the value of ficiv from the plateau at 
redshifts 3 — 5 to the local value, corresponding to an in- 
crease of a factor ~ 5. Our determinations are reported also 
in Table g] 



4 DISCUSSION 

In this paper, we have presented a new determination of 
the cosmological mass density of Civ, ilciv, at z = [1.5,4] 
based on a large sample of high resolution, high signal-to- 
noise QSO spectra which more than doubles, with respect 
to previous measurements, the covered absorption path for 
1-5 i$ z < 2.5. The main result of our calculation is that 
f2civ is no longer approximately constant in the considered 
redshift range, but shows a steady increase from z ~ 3 — 5 
to z ~ 1.5 — 2. On the other hand, it appears that the Civ 
mass density is not evolving significantly from the se redshifts 
down to the present epoch |Cooksev et alj|2009l ). 

The value of f2civ for the column density interval 12 $C 
logiV(Crv) ^ 15 can be converted into the IGM carbon 
content by mass, with the formula: 

ficiv C 



where, fib 



(5) 



CIV 

0.0224//2 2 (jPettini et al.l 120081 ) is the contri- 



7 



bution of baryons to the critical density and Civ/C is the 
fraction of C which is triply ionized, which depends on the 
assumed ionizing background. Assuming the maximum rel- 
ative abundance of Civ, that is Civ/C < 0.5, and the solar 
carbon over hydroge n abundance by mass Z c q = 0.0029 
ijAsplund et al . 2005), we obtain a lower limit at redshifts 
[2.5,4.0] of Z c > 1.4xl0" 6 = 4.9xl0" 4 Z c or [C/H] > - 
3.3. At redshift [1.5, 2.0]: Z c £ 4.6 x 10" 6 = 1.6 x 1O" 3 Z C 
or [C/H] ^ — 2.8, corresponding to an increase of a factor 
of ~ 3 with respect to the high redshift bin. 

The metallicity measured at high redshif t is in agree- 
ment with the result obtained at z = 3 by ISchave et al.l 
(2003) with the pixel optical depth method, in the case of 
over-densities above ~ 10. This is consistent with the fact 
that ficiv traces mainly the metallicity of the over-dense re- 
gions in the proximity of galaxies and it s redshift evolution i s 
linked with that of the strong systems. ISimcoe et al l |2004l ) 
found a median metallicity [C,0/H] ~ —2.8 for the inter- 
galactic gas at z ~ 2.5 with a method based on the direct 
detections of the absorption lines which takes into account 
also the upper limits. Re-normalizing our result to the solar 
abundances and the fraction of C triply ionized adopted by 
those authors (Civ/C < 0.25), we obtain [C/H] > - 3.0 
consistent with their result. 

The physical explanation for the redshift evolution of 
ficiv can be explored with cosmological simulations. We will 
devote a subsequent paper to the comparison between obser- 
vations and predictions of hydro-simulations with a detailed 
treatment of metal enrichment (Tescari et al., in prepara- 
tion). In the following, we will briefly discuss the few pre- 
dictions of f^civ present in the literature and compare them 
with our results. 

OpDcnhcimer fc Pavel l|2006l . OD06) have run cosmolog- 
ical hydro-simulations of galaxy formation meant to repro- 
duce the metallicity of the IGM, where metals are expelled 
from galaxies by momentum-driven winds (whose velocity 
is proportional to the dispersion velocity of the galaxies). 
Those winds are very effective in transporting the gas with- 
out heating it too much. Their fiducial simulation repro- 
duces consistently the star formation rate of the Universe, 
the volume averaged metallicity of the IGM and the lack 
of evolution of ficiv from z ~ 5 — > 1.5, as observed in the 
data available at the time. On the other hand, the predicted 
mass density of C increases by nearly an order of magnitude 
toward lower redshifts in the same redshift interval, due to 
the increase in the total metallicity of the gas. The constant 
behaviour of Slciv is obtained by balancing the increase in C 
abundance with galactic winds. This form of feedback heats 
the IGM causing a decrease of the C IV ionization fraction. 
In particular, the peak of the Civ ionization fraction dis- 
tribution as a function of over-density shifts toward larger 
over-densities as redshift decreases. This is consistent with 
the increase in the ratio between strong and weak C IV sys- 
tems at lower redshifts shown in Fig. [3] and the identification 
of strong systems with denser environment as deduced from 
the observed 2-point correlation function (see Section 3.3). 

To conclude, none of the models proposed by OD06 
predicts an increase of ficiv toward lower redshifts. An im- 
provement in the physics of the wind, plus the introduc- 
tion of the AGB feedback, results in a slight increase of 
Slciv in the ranges z < 0.5 and z > 5 leaving substantially 



unchanged the other values (|Oppenheimer fc Pavel liooi, 
OD08). 

An effect which has not been taken into account up to 
now is the evolution in the shape of the UV background due 
to the Hell re- ionization process expected at redshift ~ 3. 
Naively, due to the increase in the number of free hard pho- 
tons (ionizing C IV into C v) at the end of the re- ionization 
epoch, we would expect a trend opposite to what is observed: 
a decrease of the amount of C IV after z ~ 3. However, the 
detai l s of this process are still under study fe.g. iBolton et all 
120091 : iMadau fc Haardtl 120091 ) and other factors should be 
accounted for properly: the decrease in the number density 
of QSOs going towards lower redshifts and the fact that at 
low z the regions traced by the strong C IV absorbers (driv- 
ing the evolution of f^crv) could be only mildly affected by 
the cosmic UV background. 

The observed raise of the cosmic C IV mass density in 
the redshift range 1.5 — 2.5 puts a strong constraints on 
the models describing the interplay between galaxies and 
their surrounding medium, suggesting that something is still 
missing in the physical implementation of galactic feedback. 
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